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SUMMARY 


A theoretical analysis is given of the frequency and. damp ing 
characteristics of the free modes of vibrations of balanced fixed-ended 
and hinged elastic helicopter rotor blades in. hovering and vertical 
flight. Torsional vibrations, bending vibrations in flapping and in 
lagging, and coupling between the flapping and lagging motions are 
considered. Explicit methods and formulas for the calculation of 
natural frequencies and logarithmic decrements of the principal modes 
are developed, from which general conclusions are rigorously drawn. 
Flutter of helicopter blades, which may occur when the blades are 
unbalanced, is briefly considered on the basis of quasi-stationary flow, 
and staple criterians are derived for the stability of the coupled 
torsional and flapping vibrations in such cases. 


INTRODUCTION 


The purpose of this investigation is to present under one cover, 
in a brief and simple fashion, a comprehensive analysis of the frequency 
and damping characteristics of elastic helicopter rotor blades in 
hovering and vertical flight when they perform small vibrations about a 
state of static equilibrium. Effects which appear not to have been 
thoroughly investigated heretofore are considered in detail here. These 
include: Consideration of the effect of boundary conditions on the 

centrifugal contributions to the natural vibration frequencies of 
rotating blades in bending; the effect of free torsional vibrations on 
the flapping and lagging vibrations of mass-balanced blades; the effect 
of Coriolis, centrifugal, and aerodynamic coupling between the flapping 
and lagging vibrations of a helicopter blade; and a comparison between 
aerodynamic and internal damping in the principal modes of rotating 
beams in bending and in torsion. 
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The analysis and general formulas have been developed for blades 
having any taper and cross-sectional distribution, with the restric- 
tion that the shear center and the center of gravity of a cross section 
coincide (as is usually the case, for example, for a single— tubular- 
spar section). Moreover, the case of mass-balanced blades, that is, 
blades whose cross-sectional centers of gravity coincide with the aero- 
dynamic centers, has been emphasized, since attempts are usually made 
in practice to achieve such a condition. (See reference 1. ) Neverthe- 
less, the vibrations of unbalanced blades, characterized by coupling 
between flapping and torsion with the resulting possibility of flutter, 
have been briefly analyzed. 

In the case of hinged blades, it has been assumed that the flapping 
(horizontal) and lagging (vertical) hinge axes are intersecting and 
perpendicular to the blade axis, so that no change in pitch angle of-a 
blade is caused by either flapping or lagging. 

The general procedure in this analysis consists first in setting 
up the equations of motion of a rotating helicopter rotor blade in 
hovering and vertical flight in torsion and in bending in two mutually 
perpendicular directions (flapping and lagging) and then in solving 
these equations either exactly or else approximately by the Rayleigh— 
Ritz method. 

Unless otherwise stated, the details of the mathematical deriva- 
tions have been almost entirely omitted here for the sake of brevity; 
however, these details have been worked out by the present author in an 
unpublished paper entitled "A Theoretical Analysis of the Elastic 
Vibrations of Fixed— Ended Helicopter Blades in Flight." 

This work was conducted at the Polytechnic Institute of Brooklyn 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The author wishes hereby to express 
his appreciation for the kind aid given him by Dr. Hans J. Reissner. 

The author also expresses his gratitude to Messrs. Charles Chin, 
Constantine Dolkas, Francis G. Hinchey, and Frank J. Romano for their 
assistance in the numerical calculations. 


SYMBOLS 

A cross-sectional area of a blade 
Aq - value of A at root of a blade 
B downwash factor, equation (2b) 
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dimensionless parameter 



length, of chord of a blade section 

value of c at root of blade 

value of profile— drag coefficient of blade section when 
absolute angle of attach is equal to 0 O + 0p 

value of profile-drag coefficient of blade section at 
zero absolute angle of attack 

dimensionless distance between midpoint and shear center 
(here, also center of gravity) of a blade section 
(fig- 1) 

modulus of structural blade material 

distance between point of attachment of blades and rotor 
axis of rotation; eccentricity- 

shear modulus of structural blade material 

acceleration due to gravity 

. internal damping coefficients in flapping, lagging, and 
torsion, respectively 

line of intersection of plane (x,y) of rotation and plane 
of cross section of blade (fig. 1) 

structural moments of inertia of a blade section for 

bending in flapping and in lagging planes, respectively 

values of Ip and I 2 at root of a blade 

mass moments of inertia per unit of blade length of a 
cross section about principal axes through center of 
gravity (here, also shear center) 

mass polar moment of inertia per unit of blade length of 
a blade section about principal axes through center 
of gravity (here, also shear center) 

values of 1^ 


and J at root of a blade 
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i - \T=i 

j 


K i 

*2 

k 


torsional moment of resistance of a "blade section 

dimensionless distance between center of gravity (here, 
also shear center) and aerodynamic center of a blade 
section (fig. 1) 

dimensionless bending-stiffness parameter ( EI lo/ aA o^ 2 
dimensionless bending-stiffness parameter (^EIpo 

ka 2 


constant in the relation c 


do 


( Oa °)o + 


m 

Mi 

Mg 
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length of blade 

number of blades in rotor system 
dimensionless torsional-rigidity parameter 

_ ~ ^min)o 


(“0/ n %o l2 ) 


I P° 


subscript indicating a given principal mode of vibration 

P* c o 3 / 


dimensionless parameter jt 


"^PC 




complex frequency; if p = — R ± ico (R and a> real) 
then co/2ir is the natural frequency in cycles per 
second, while 2k & is the logarithmic decrement 


Q 

q. 

R 
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dimensionless parameter 


pl 2 c 0 2 


■‘pc 


dimensionless complex frequency (p/ft) 
tip radius of a blade from axis of rotation 

t 

local radius of a blade element from axis of rotation 
time 
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v 


w 


dimensionless "bending (lagging) deflection in the direction 
of a principal axis y^_ (fig* 1) of a blade section 

gross weight of helicopter 


v 


v(t) 


dimensionless bending (flapping) deflection - in the direction 
of a principal axis z^’ (fig. 1) of a blade section 
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0(t) = 

0 0 
e l 

1 

P 

cr 


induced downwash velocity in steady state 


local slope of a deflected blade in flapping in the 
steady (static) state 


local absolute an^ie of attack of a blade section 

aerodynamic and internal (structural) logarithmic decre- 
ments , respectively 

eccentricity ratio (e/z) 

local twisting angle of a blade section 


§z zz 0 ( 1 ^) 


■'O 


design, or initial, blade a n gle (fig. 1) 

angle between zero— lift line and principal axis y^_ of 
a blade section (fig. 1) 

dimensionless distance along a blade, measured from root 
air density 

average density of blade material 
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t(§) 


dimensionless centrifugal— force parameter 



A_ r 
A o 2 



n angular speed of rotor system, radians per second 

(jo natural frequency of Titration of an elastic rotating 

tlade in a given principal mode 

co cn value of m if tlade had no tending stiffness, that is, 

if Z 1 = = 0 ■ 

oi en value of % if Cl = 0 

t = <L 

• = £L 
St 


* 


BASIC EQUATIONS 


The tasic differential equations for the small flapping, lagging, 
and torsional vibrations, respectively, of a helicopter tlade in 
hovering and vertical flight about its position of static equilibrium 
can be written in the following nondimens ional form: 
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( 1 c) 


Equations (la) and (lb) represent the equilibrium of the elastic, 
internal damping, centrifugal, inertia, and aerodynamic loads at any 
point | per unit length along a blade. In equation (la) these loads, 
as well as the dimensionless displacement w, are in the direction of 
the principal axis z-j_ of a cross section of a blade (see fig. 1 ). 
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In equation (lb) the loads and the displacement y are in the direction 
of the other principal axis of a cross section (fig. 1). If the 

blades are hinged then w and v represent the sum of_the elastic 
deflections and the components of the rigid-body displacements in the 
directions of the principal axes z-j_ and yj_, respectively. 

The aerodynamic loads here have been derived, as in reference 2, 
by use of the three-dimensional Kutta— Joukowskl theorem for quasi— 
stationary flow. The induced downwash, however, has now been assumed 
as constant throughout the rotor disk and as derivable by the simple 
momentum theory. (See reference 3.) The only difference, in fact, 
between the hovering state and the state of uniform vertical climbing or 
descending In equations (la) to (lc) lies in the values of the downwash 
velooity w-j^q and of the downwash parameter B. 


By the momentum theory, the values of and B for a heli- 
copter climbing vertically at a constant speed of v c can be shown to 
be: 



y 


where the drag of the helicopter in climbing has been assumed to 
be kpW, and V is the gross weight of the helicopter. 

It has been assumed, in the derivation of equations (la) to (lc), 
that the profile-drag coefficient c^ of a blade section varies 

parabolically with the angle of attack. 


From the condition that the static thrust must support the gross 
weight of the helicopter, the required value of the absolute pitch 
angle (0 O + of a blade, assumed constant- along its length, is 
found to be: 


9l = 


w • 

jtpmc fl 2 R 2 1 
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The mechanical significance of each of the terms appearing In 
equations (la), (lb), and (lc) is as follows. The terms proportional 
to Eq, Eg, ° r Mi represent the elastic resistance, while those 
proportional to g w , g T , or Sq represent internal damping. (See 
reference 4.) The terms ( tv 1 )* and ■£- 0 o v in equation (la), the 

terms ( tt 1 ) 1 , — v, and — 0 o w in equation (lb), and the Mo— term in 
Ad Ad 

equation (lc) are due to the centrifugal loads. The terms with 


v. 


and 0 represent inertia loads, while the terms 

t 


and 


2 A. Vr 


w 


p A w Q « T 
A 0 i a 


In equations (la) and (lb) 


-2— — give the Coriolis effect. 

Aq l a 

the terms proportional to C represent the aerodynamic loads, while in 
equation (lc) these loads are represented by the terms proportional 
to Q and to P. For mass-balanced blades, Q = 0. The terms in w 

-4- -A- e 


and 9 give the effect of induced downwash. Finally the term 
in equation (lb) represents the effect of the weight of a blade. 


n 5 * Ad 


TORSIONAL VIBRATIONS 


From equation (lc) it can be seen that if the blades are mass- 
balanced, implying Q = 0, then the torsional vibrations will not be 
coupled to the flapping vibrations. This case of balanced blades is the 
one which is treated in the present section on torsional vibrations and 
in the following two sections on vibrations in flapping and lagging. 

Frequency characteristics.— By putting 


Hi,*) - 8(t)» pt 


W 


into equation (lc) with Q = 0, and neglecting damping for the present, 
the following equation is obtained for 0(|) and q = ^: 
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0 = 0 


(5) 
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If the "blades have a uniform, cross-sectional distribution and have 
fixed pitch at the root (0=0 at | = 0), then the following solution 
is obtained for the principal torsional modes of vibration: 


0 n (|) = sin n J | 


n odd, n ^ 1 


( 6 ) 


V - M 1 t + «2 ~ % f n2 


The term 0 n ( £ ) represents the principal modes of deflection, 
while \/-q 2 is the natural frequency in cycles per second in 

"n+i rt " 

the -^=th mode, if 0 is in radians per second. 

It is significant to note from the results of equations (6) that, 
since M2, which represents the effect of the centrifugal torque, is in 
actual cases negligible, the natural frequencies, as well as the modes 
of vibration, of a rotating helicopter blade in torsion are virtually 
the same as those of the same blade when it is stationary (fl =0). 

This conclusion is valid for tapered, as well as for rectangular, blades 
(cf. appendix A). 

If the blades have a variable cross-sectional distribution, then 
the natural frequencies of the various modes of vibration can be deter- 
mined by the procedure, based on the Rayleighr-Ritz method, shown 
explicitly in appendix A. 

Damping characteristics.— Let the natural frequency of the "kth" 

principal mode of vibration without dancing be o^, and let the corre- 
sponding mode of deflection be 0 no (|). Then it” may be assumed that for_ 
this mode with dancing the deflection shape will be given approximately 
by 0 no ( | ) - Therefore, with the assumption 


SntS-t) - 1 ’W^ ePt 


for the damped vibrations of the kth mode, where p may now be complex, 
differentiation of the integral in condition (Al) of appendix A with 
respect to b leads to a quadratic equation in q whose roots (with M2 
neglected) indicate the following values of the natural frequency 0^4, 
aerodynamic logarithmic decrement 5^, and internal logarithmic 
decrement 5 in of- the kth mode: 
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where 







ffc) 3 ^ g ft> ^ 

?- <W^ ) « 

T>0 


Be ay 
2 


no 

n 


( 7 ) 

( 8 ) 


( 9 ) 


(7a) 


Equation (7) gives the correction, due to damping, in the natural 
frequencies. This correction will usually he found to he negligible in 
all inodes above the fundamental. 

Equation (8) shows that the aerodynamic logarithmic decr emen t in 
any given principal mode has an order of magnitude of — and varies 

approximately inversely as the natural frequency. Thus, the aerody nami c 
logarithmic decrement decreases with the mode of vibration. Equation (9) 
shows that the internal logarithmic decrement is, on the other hand, 
independent of the natural frequency and therefore remains the same in 
all modes. This result for the internal damping is due to the negligible 
effect of the centrifugal twisting couple in any cross section of a 
rotating helicopter blade. 

Since in general g Q « P, it follows from equations (8) and (9) 
that in the fundamental, and possibly second, mode the aerod ynami c 
damping will be greater than the internal damping, in the high modes, 
however, the internal damping will predominate. 
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It should he noted that, although equation (9) indicates the same 
internal logarithmic decrement for stationary as for rotating beams, 
equation (8) shows that there will be no aerodynamic damping if the beam 
is stationary (5^ = 0 if fl = 0). 

Numerical exangle .— Consider a blade of constant cross section along 

the span, with the data of figure 2. Then e_.(£ ) = sin n £ | for the 

G+l 2 

-g— ‘ th mode. Moreover, it is found that Mj_ = 37. 0 and P = 7.35. In 

virtually all cases. Mg = 1. Here let =_0.7. Finally let g fl = O-O5. 
From equations (6) it is found that (damping neglected) ' 

pp » \j 91. On 2 + 0.7 a \/ 91. On 2 = 9.5811 

The fundamental (n = 1) natural frequency in torsion is thus quite 
high, being here almost 10 times the rotor angular speed. From equa- 
tion (7a), Tj nn = 7*35 ^P* 2 5 + — + 0.239n. For the fundamental 

mode (n = 1), equation (7) therefore gives: 


^p = ^91.0 - (2.82) 2 = 9.10 


indicating thus a decrease, due to damping, of about 5 percent in the 
fundamental frequency. The effect of the damping on the natural 
frequencies in the higher modes (n = 3, 5, . . . ) will be much smaller, 

as can be seen from equation (7) with the expressions for and ^p 

here in terms of n. From equation (8), it is found that = I.78O 

1.20 

and 5^ % ■ " n ■■ for n ^ 3. Moreover, 6^ = 0.1570 for any mode. 

Thus in the fundamental mode the aerodynamic damping is here over 
10 times as great as the internal damping. In the second mode (n = 3), 
however, 6 a is reduced to less than one-fourth of its value in the 
fundamental mode, while 5^ remains the same. In the fifth higher 
modes (n >9)* the internal damping exceeds the aerodynamic damping. 


BENDING YIBRATIQNS IN FLAPPING 


In considering bending vibrations in flapping the small coupling 
between the flapping and lagging vibrations will be neglected. Moreover, 
it will be assumed here that the blades are mass-balanced. In t ha t case 
the torsional vibrations can be determined first, as In the preceding 
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section* "by obtaining 9 explicitly as a function of f and t. This 
known function can then be substituted for 9 into equation (la)* which 
is the basic equation for the motion of 3 blade in flapp ing . 

Frequency characteristics .— The natural frequencies of the free 

vibration modes of a rotating blade in flapping are determined by the 
equation 


El & w ”J' (TV,,,+ t^ = 0 Cla,) 

The effect of dairying is not included in equation (la» ). By putting 

w(fijt) = v(g)e pt 



equation (la*) becomes: 

K i(^~ ^ “ C™*)’ + 4 2 w = 0 (10) 

The values of q 2 determined by equation (10) for any principal 
mode of vibration can be expressed approximately in the form (see 
reference 5) : 

2 2 2 

Sl = 9 en + 4 on (“> 

O p 

where is the value of q n when the beam is stationary (0 = 0)* 

and q cn 2 is the value of q n 2 when the blade is hypothetically 
rotating without bending stiffness (fo = 0). The values of a 

en 

represent the effect of the elastic resistance of the blade* while the 
values of q^ represent the effect of the centrifugal loads. 

The following remarks conceming relation (11) are pertinent. It 
can be shown (appendix B) that relation (11) would be exact if and only 
if the centrifugal loads would have no effect on the modes of deflection 
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of a rotating "beam. It can also be shown (appendix B), however, by a 
single consideration of boundary conditions, that equation (11) cannot 
be an exact general statement f<jr either- fixed— ended or hinged blades 
and that the centrifugal loads must therefore in general have some 
effect on the modes of deflection. It should be further noted that the 
values of q^ must, according to relation (11), be independent of the 
boundary conditions of a blade (cf. appendix in reference 5)- It can 
be shown (appendix B), however, that there are theoretically no exact 
values of q^ 2 for a blade fixed at the root. (This may be considered 
as a simple proof that relation (11) cannot be exact for fixed-ended 
blades. ) 

In spite of the fact that equation (11) cannot be an exact general 
statement, this relation gives sufficiently accurate results for 
practical purposes. This is due to the fact that^ although the centrif“ 
ugal loads do have an effect on the modes of deflection of a rotating 
beam, this effect is ordinarily small both for fixed-ended and hinged 
blades. (See references 6 and 7- ) As a numerical check on the accuracy 
of relation (11), the values of the natural frequencies of a hinged 
blade were calculated by solving equation (10) without the use of 
relation (11). Approximate solutions for the lower modes (see appendix B 
for details of the method) showed, that- the natural frequencies as func- 
tions of the parameter could to a high degree of approximation be 

represented in the form of equation (ll). Moreover, a numerical check 
(appendix B) was made on the effect of the root conditions of a blade 
on the centrifugal contributions q^ 2 to the negative squared natural 
frequencies, and this effect was found to be practically negligible. 

The values of q^ 2 in relation (11) depend only on the mass 
distribution — (| ) of a blade along its length and are the same 

( aA )o ? 

for a fixed-ended as for a hinged blade. The values of q^ , on the 

other hand, are in general, directly proportional to the parameter Kp 
of a blade and depend on the boundary conditions at the blade root, as 
well as on the cross-6 ectional distribution (i.e., plan form). 
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For "blades of constant cross section, the Yalu.es of q^ 2 and q^ 2 
in the various principal modes are as follows: 

■ fixed— ended "blades: 

= 22 - 3 ^> 481, . . ., n)^" (approx. ), . . . 

n odd, n ^ 5 


— q^.^ 2 = 1, 6, . . . , n(2n 1), . - . 

n a positive integer 

hinged blades: 

-<l en 2 /% = °» 2 3 6 j • • *.» *Vn + 0*” (approx. ), . . . 
n a positive integer 

— q^,^ = lj 8, . . », n(2n — 1), • > • 
n a positive integer 

For blades of variable cross section a simple procedure, based on 
the Rayleigh-Bit z method, of calculating the frequency characteristics 
is given explicitly in appendix B. 

The values of the frequency ratios co/ft in the various inodes of 
vibration of a blade in flapping are evidently functions of only one 
parameter, K-^, in addition to being .dependent on the plan form of the 
blade. Typical values of Kq are usually small,. For the single- 
tubular-spar section in figure 2, for example, E-^ = 0.00400. Conse- 
quently, the actual contribution of the elastic resistance of a heli- 
copter blade to the natural flapping frequency in the fundamental mode 
is small in comparison with that of the centrifugal forces. In the 
higher modes, however, the relative importance of the elastic resistance 
rapidly increases. The natural frequency of a fixed-ended blade in any 
principal mode will, of course, be higher than that of the same blade 
when it is hinged. If Ej_ is of the order of magnitude of unity this 
difference will be especially great in the lower modes. However, 
when EjL « 1 (as in the' tubular-spar section, fig. 2), then the 
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natural frequencies of a fixed— ended Hade are only slightly higher than 
those of a hinged "blade. These conclusions are, in part, illustrated 
"by figure 3 and figure k, where numerical results are plotted for a 
uniform "blade of low "bending stiffness hinged at the root and fixed at 
the root, respectively. 


Rigid-body oscillation as a mode of vibration.- If the blades are 


hinged at their roots, and if their flapping hinges pass through the 
axis of rotation of the rotor system, so that ^ = g, then the deflec- 


tion shape w = ag, where a is a constant, will be an exact solution 
of equatiop. (10) with q 2 = —1, regardless of the plan form of the 
blades. This shows that, if damping and small coupling effects are 
neglected, then a rigid— body oscillation of a blade in flapping will be 
an exact fundamental mode of vibration. In this mode, co = fi exactly. 
If, however, aerodynamic damping is considered, implying addition of 

X 1 C / w \ 

— — Bw ! to the left side of equation (la'), then it 

At / 


the term 


° T c^Vi 


can be seen that equation (10) will no longer be exactly satisfied 
by w = ag. In this case, the rigid-body oscillations will be an 
approximate, but not exact, mode of vibration of the blades. 


Damping characteristics .— Putting w(g,t) = w(g)e£’ fc Into equa- 
tion (la») with the term C ^ - Bw) added to the left side, the 

following equation for w(g) and q Is obtained: 



(tw 1 ) 1 + — q 2 w + C — — q (vr — Bwj^ = 0 

A o Z C ri ' ' 


By mathematically transforming this equation into"a stationary condi- 
tion, 1 and applying the Rayleigh method with the assumption that the 
damping has a negligible effect on the principal modes of deflection, 
the following relations are obtained for the effect of the aerodynamic 
damping on the natural flapping frequency of any mode (denoted by 
subscript n) and for the aerodynamic logarithmic decrement- 5 an in 
any mode In flapping: 2 


1 It is interesting to note that it is here possible to derive 
rigorously the stationary condition of an integral for a nonconseryatlve 
system after. making the substitution w(g,t) = w(g JeP^. 

O 

In terms of the value of the complex frequency q , for any mode 
with damping, equation (13) Inplies; ntt 

«n 4 " “(I WV) ± 1 TT 
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“>nd = 

a 



( 12 ) 


5 


an 


ynn/ g hn 

%d/n 


(13) 


where tOp^ is the natural frequency for the damped mode, 
corresponding natural frequency without damping, and and 

are constants for each mode, defined as: 


is the 
c Wl 


7 = C 

'nn 1 


-2_ H w *(t) it 

c 0 j no a-5 





a = 
nn 





Here, w no (|) i s the mode of deflection without dancing. 

The Internal logarithmic decrement, derived and discussed in detail 
■in reference 4, is : 


^in ~ ag w 



( 1*0 


From equation (12), which gives the effect of the .aerodynamic 
damping an the values of the natural frequencies, it follows that this 
effect decreases with the mode and will in actual cases be negligible 
in all principal modes with the possible exception of the fundamental. 

In the fun damen tal mode, in fact, the damping may in some cases (namely, 
those cases for which the parameter C is sufficiently large) be so 
heavy that this mode will consist of an unoscillating decaying motion 
instead of an oscillating motion. 
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Since the actual value of -l' nn / a < an will he roughly the same (of 
the order of magnitude of C/2) for any mode., it follows from equa- 
tion (13) that 8 an will decrease with the mode, in virtually inverse 
proportion to the natural frequency. Comparison of actual orders of_ 
magnitude (0(C) = 1, 0(g w ) = 0.05) shows, according to equations (13) 
and (14), that except in the very high mode?, where + 5^ £ 5^ £ jrg v , 

the aerodynamic damping will exceed the internal damping. The latter is, 
in fact, negligible in comparison with the aerodynamic damping in the 
fundamental mode. Because of the aerodynamic effect, the flapping 
vibrations in the lower modes will in general he highly damped. 

Numerical example .— For a uniform blade with the data of figure 2, 
one finds: K-j_ = 0.00*1-00, C = 1.7^, and B = 0,532 (equations (2a) 
and (2b) with v c = 0 and R = 2). For fixed— ended blades, it is 
permissible to substitute w n0 (g) = 5^ — + 6S; 2 into equations (12) 

and (13) for any mode. 3 For hinged blades, the substitution w nQ ( £ ) = £ 

may be made in these equations. It is thus found ^assuming y = 

that ^22: - = 1.105 when the blade is fixed-ended, and = 0.995 

Tin 2.31 o-nn 

when the blade is hinged. Moreover, for fixed-ended blades, 

^2. = \/ 12.36 x 0.00*4- + 1 = 1 . 02 k- ; for hinged blades, 

= \j 0 x 0 . 00 k + 1 = 1 in the fundamental mode without damping. 

With aerodynamic danping taken into account, these fundamental frequen- 
cies change thus: 

= \f( 1.024) 2 - = 0.862 

\ sz / fixed-ended v \ 2 / 

- *•** 


Thus aerodynamic danping diminishes the fundamental frequency here by 
about 13 percent for both the fixed-ended and the hinged blade. 

In the second and higher modes, however, the effect of the damping 
on the natural frequencies is here negligible, since the value 
of rapidly increases (over 6 in the second mode), while the 

^Since the value 7 Taj J^ xai ioes noi: Yar ¥ greatly with the func- 
tion w nQ (|), it -follows that, in actual calculations, rough assumptions 
for w no ( | ) in equations (12) and (13) will suffice to give sufficiently 
accurate results for practical purposes. 
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value of f 7 nn / 2 a im 'j 2 remains approximately the same, at out l/4 here. 
The total logarlthml o decrements (8^ + 8^) for "both the hinged and 
fixed— ended blades, with g v assumed as 0 . 05 , are shown in figures 5 
and 6, respectively, plotted against the principal mode of vibration. 

For comparison with internal damping, the aerodynamic logarithmic 
decrements are also plotted separately in these figures. 

Effect of torsional vibrations on flapping vibrations .— The free 
torsional vibrations of a mass-balanced blade in any principal mode 
have an effect on the free flapping vibrations of the blade similar to 
that which an external damped harmonic load would have, if the frequency 
and damping of this load are the same as the frequency and damping of 
the corresponding torsional mode. This effect can be determined by 

p t 

substituting the known function 9 = 0 n ( £)e 11 for any principal mode 
in torsion (cf. TORSIONAL VIBRATIONS) into equation (la) and then 
solving for w(|,t). By this procedure it can be shown that corre- 
sponding to each principal mode of vibration in torsion of frequency cj^ 
and logarithmic decrement S n there will appear in general a component 
of vibration in flapping of the same frequency and logarithm! c decrement. 
It can be shown, however, that in all actual cases the ratio of the 
amplitudes of these " quasi— forced" flapping vibrations to those of the 
corresponding "quasi— forcing" torsional vibrations will be of the order 

of magnitude of - — ^ + — — — i — . Since in practice ^ » 1 

K/n) 2 it i <%/ n n 

(cf. TORSIONAL VIBRATIONS) it follows. that the relative amplitudes of 
the quasi— forced flapping vibrations will be so small that the effect 
of free torsional vibrations on the flapping vibrations of mass— balanced 
blades may In practice be neglected. 


BENDING- VIBRATIONS IN LAGGING 


The characteristics of the lagging vibrations of a helicopter 
blade, virtually in the plane of rotation, are discussed in this section. 
Coupling effects are here neglected and are discussed in the next 
section. 

Frequency characteristics .— According to equation (lb) the natural 

frequencies of the free uncoupled lagging vibrations of a blade, with 
damping at first neglected, can be determined from the following 
ordinary differential equation for v(£) and q: 



(tV* )* + -^(q 2 - l)v = 0 


(15) 
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A comparison of equation ( 15 ) with equation (10) shove that 
I2 • Ii 

if — — (|) = — — (|), as will ordinarily he the case, then the following 
-*■20 x lo 

simple relation exists hetween the natural frequencies in lagging and 
flapping: 


oij 2 = cOf 2 (K2) - n 2 


(l6a) 


where is the natural frequency (in cps) in lagging for any principal 

mode, 05 f (K2) is the natural frequency of the corresponding mode in 
flapping with replaced hy Eg, and Cl is the rotor angular speed 

(in rps here). Relation (l 6 a) is valid for fixed-ended blades as well 
as for blades with both flapping and lagging axes hinged. 

From equation (11) for flapping, it is evident that relation ( 16 a) 
can also be expressed in the following form for any principal mode 
(characterized by subscript n)': 


“^nZ “enZ 2 + “cnZ^ 


where 


2 2^-2 
“enZ = <°enf f 


(D. 


cnZ 


= o>. 


cnf 



(16b) 

(16c) 

(l6d) 


The subscripts Z and f refer, respectively, to the values in laggixg 
and in flapping. Thus the elastic contribution to the square of the 
lagging frequency of- any principal mode is the same as the elastic 
contribution- to the square of the corresponding flapping frequency, 
except that the former is proportional to Kq while the latter is 
proportional to Kj_. This fact is represented by figure 7* From 
equation (l 6 d) it follows that the centrifugal contribution to the 
square of the lagging frequency of any mode is less by a fixed 
amount (fi 2 ) than the centrifugal contribution to the square of the 
corresponding flapping frequency. 

From equation ( 16 a) it follows that if Eg = (as in single- 
tubular-spar sections) then the lagging frequencies of all principal 
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modes will "be lower than the corresponding flapping frequencies. In 
the fundamental mode* where r; ft (because of the actual low values 
of the bending— stiffness parameter Ei), the lagging frequency will be 
very low. In the higher modes , however , the difference between the 
flapping and lagging frequencies , if E -]_ = becomes relatively less 

and less. If, on the other hand. Eg » Kp (as for structural airfoil 
sections) then the natural lagging frequencies will be greater t han the 
corresponding flapping frequencies, often even in the fun damen tal mode. 
In that case, in fact, the following simple relation will be valid for 


o So o 
the high modes: (a-f- ts — . 


1' 


It may be remarked that if = Eg and 


12 Ii 

r(5 )s r 

x 2o Ilo 


(l), then 


the principal modes of deflection of the free und amp ed uncoupled lagging 
vibrations will be exactly the same as those of the corresponding 
flapping vibrations, provided of course that the root conditions of the 
blades be the same (e.g., both fixed— ended or both hinged) in both 
flapping and lagging. ^ 


Effect of a small eccentricity .— If the blades are hinged at the 
rotor axis of rotation then in the fundamental flapping mode, = ft. 

In this case, equation (l6a) inplies a>j = 0, signifying the absence of 
a centrifugal restoring force. Such a situation Is easily remedied by 
attaching the blades at a small distance (eccentricity) e away from, 
instead of at, the axis of rotation. The effect of such an eccentricity, 
for both hinged and fixed— ended blades, can be expressed by the relation; 


2 2 w 2 

“ = “c=0 + E !» 


(17) 


vliere cd and ft are In cycles per second, e - — , 

& 


E H I X»2( S) [ 
0 




a s 


^2 2 U) cLS 


^The condition E]_ = E2 is 
effect of the centrifugal loads 


necessary only because of the small 
on the modes of deflection. 
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and X(| ) is the mode of- deflection in flapping. Equation (17) is 
valid for any principal mode in either flapping or. lagging. 


In accordance with the Rayleigh method, the natural frequency 
in flapping of a flexurally weak cable (Kp = 0) with e = 0 

can be expressed approximately 5 in the form: 




where 


t(|) 



Hence equation (17) can be written approximately in the form: 



The coefficient of (^ c f)e=o 2 in equation (18) is almost independent 
of the mode of deflection X(|) and is greater than unity for any mode. 

It is evident from equation (18) that an eccentricity of attachment 
of the blades increases the natural frequency of any mode in either 
flapping or lagging and that this increase is .greater the higher the 
mode, since <a C f 2 increases with the mode. However, for the reason 
already explained, this eccentricity effect is especially important in 
the fundamental lagging mode of a hinged blade or of a fixed-ended 
blade with a low bending stiffness in the plane of rotation. 


For the fundamental mode, it is permissible to substitute— X(| ) = 5 
for hinged blades or X(l) = — 4^ + 65 for fixed-ended blades into 

^The approximation in the equation consists in neglecting the 
effect of the centrifugal loads on the modes of deflection. 
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equation ( 17 ) or (l8). 

IT 

is thus found that — 

Gr 

fixed— ended "blade. 

Dancing characteristics.— In a manner analogous to the derivation 

of equations (12), ( 13 ), and (14) for flapping, the following values 
can "be obtained for the natural frequency and for the aerodynamic and 
the inte rna] logarithmic decrements in any damped principal lagging 
mode (characterized "by subscript n): 



If — (I) = 1, then for equation (17) it 
A© 

= 1.5 for a hinged blade and ^ = I.558 for a 



and Y n (|) is the undamped mode of deflection in lagging. The value 
of the ratio 7 nri /% n is almost independent of the deflection 
shape Y n (|). 

Since in actual cases the quantity (j^/2 a ^ 2 will be negligible 

in comparison with (%. 0 /^^ ^ follows from equation (19) that the 

effect of damping on the lagging frequencies is in general negligible 
In al i modes. The logarithmic decrements in lagg in g will vary with the 



2k 


NACA TN 1999 


mode of vibration in a manner quite similar to t hairpin flapping (see 
BENDING VIBRATIONS IN FLAPPING) . It must - be observed, however, that the 
magnitude of the total (i.e. , aerodynamic plus internal) damping is, in 
the lower modes, much smaller than that in flapping. The value of 5^ 

( c do)l U^nf 

in lagging is, in fact, of the order of — times that in 

it <%2 

flapping. In the fundamental mode, this ratio may be about - l/30. The 
internal damping in 'lagging is at least of the same order of magnitude 
as the aerodynamic damping in the fundamental mode of fixed— ended blades 
and is therefore relatively much more important in lagging than in 
flapping. 

Numerical example .— For the uniform blade with the data in figure 2 
and e = 0.05, the values of the undamped natural frequencies in lagging 
are plotted against the mode of vibration in figures 8 and 9 for hinged 
and fixed-ended blades, respectively. The logarithmic decrements in 
lagging with (Odo^ = °*°3 31111 8v = °*°5 ar 0 plotted in figures 10 

and 11 against the modes of vibration. 

Effect of torsional vibrations on lagging vibrations .— The free 
torsional vibrations of mass-balanced blades have an effect on the free 
lagging vibrations which is quite similar to thair which they have on the 
free flapping vibrations (see BENDING VIBRATIONS IN FLAPPING). This 
effectj which is mathematically represented by the terms with Q in 
equation (lb), can, however, be shown to be negligible, because of the 
very low order of magnitude of the relative amplitudes of the lagging 
vibrations induced by the torsional vibrations. 


CENTRIFUGAL, CORIOLIS, AND AERODYNAMIC COUPLING 
BETWEEN FLAPPING AND LAGGING 


In general, there will be centrifugal, Coriolis, and aerodynamic 
forces which couple the flapping and lagging vibrations of a helicopter 
blade. These forces are represented mathematically in equations (la) 
and (lb) as follows: The terms 0 o v and — 9 0 w represent the 

centrifugal coupling loads per unit length of the blade; the 


terms 2 


A w^ 
^ l 


— and 


fl 


finally the terms — C 


T#o 




V 

— represent - the Coriolis load; 


+ G • 




LLo 


nr 


I 

n 
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and C y y-^20 o + Sq — 2 — B f v ) re P reseirfc fke aerodynamic 

coupling loads. ^ All these loads are evidently of a relatively low 


order of magnitude "because of the actual low values of Q 0 and 
w 0 » 

of — — . Consequently., for a first approximation, these coupling terms 


may, as in the preceding two sections, he neglected in determining the 
frequency and damping of any principal mode in either flapping or 
lagging. The effects of the coupling on the frequency and damping 
characteristics of. any mode can then "be taken into account "by making 
the corrections described in the following discussion. 


Effects of coupling on frequency and damping characteristics .— By 
putting 


v(S,t) =‘w(^® pt 

= v(l)e pt 


into equations (la) and (lb) it is possible to derive the following 
expressions for the complex frequencies and for any 

principal mode (characterized by subscript n) in flapping and in 
lagging, respectively, with the simultaneous coupling loads taken into 
account: 


Wl 


L nfo 


2 ^nf o ^fen ^n^^^n ^n^ g ^n^n^plnfo 

- Vlo 2 )( aa hVo ♦■’’n) 


Sa 2 »o 2 


(22a) 


^nZl SiZc 


2 SiZo (^n %^ 2 ^n 2 ^n^o^n^-nZo ^n Q c 

^ a nSiZo^nZo — ^nfo^^n + a n(^nZo + ^nfo)J 


(22b) 

b It may be noted from equations (la) and (lb) that the aerodynamic - 
loads directly oppose the Coriolis loads but do not. in general cancel the 
latter. 
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where q^Q and q nlo are the cougjlex values of q in the corre- 
sponding uncoupled , though damped, modes in flapping and lagging, 
respectively, and where cc^, c^, . . . are real constants for each 
uncoupled mode defined thus: 


a = 
n 


W 2 dt 
Aq u 


a n s 


t Y » 11 


Pn s 


— W Y d| 
Aq n v n 


K = 


n 


£!sl wt d4 

Aq l " nVn ' 


d 0 


± m0 f ^( 28 ° + * 01 " 2 




7n 


■° I fl V n 2 ^ 

o 


Here W n (£) and Y n (|) are the modes of deflection of the uncoupled 
flapping and lagging vibrations, respectively. 


An advantageous property of equations (22a) and (22b) is that, 
according to the forms of these equations, the numerical-values of the 
corrections to q n f Q and <3. n ^ 0 will be rather insensitive to the 
forms W n (|) and Y n (|) of the deflection shapes substituted into 
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these equations* For purposes of a quick calculation it is therefore 
suggested that the substitutions V n (£ ) = Y n (£ ) = | 

and ¥ n (|) = Y — 4|3 + 6g 2 be made for any mode for hinged 
blades and for fixed-ended blades, respectively. 

From equations (22a) and (22b) it is found that in general the 
order of magnitude of the corrections, due to coupling, in the natural 
frequency and logarithmic decrement of any principal mod^ in flapping 

or lagging will be that of second powers of 0 O and - 2 — . Thus, the 

coupling will have what may be considered a second— order effect on the 
frequency and damping characteristics of any mode. This effect may, 
nevertheless, be quite appreciable in the fundamental lagging mode, 
since the damping decrement (i.e., the negative real part of q) in 
this mode is generally very small, while, if the bending stiffness 
(represented by the dimensionless parameter Kg) of the blade in lagging 
is small (as in the tubular— spar section, fig. 2) the natural frequency 
in this mode will also be low. 


Numerical example .—. Consider a uniform blade with the data in 
figure 2. Also let Sp = 0.02, e = 0.05, and g T = g w = O.O5. Then 
from equation (3), 0 O = 0.1273* 

(a) Fixed-ended blades; For the fundamental mode in flapping (see 
Numerical example under BENDING- V1BBATI0NS IN FLAPPING and equation (17) ), 


+ i M 1 + 22.36 X 0.004 + 0.05 X 1.558) - 


^ 1 . 105 \ 2 


= -0*553 + 0.907i 

For the fundamental mode in lagging (see figs. 9 and 11), 


0.178 X 0.357 . 

*110 = £ + °' 35Ti 


-0.0101 + 0. 357i 


Putting W^g) = V x (|) = ^ - 4g 3 + 6g 2 and 
that: a = = P 1 = 2.31, = 0.741, and 


Y = | + e, it is found 
7j_ = 3.43. 
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W n * Jt 7 — 

Assuming — — (I ) = 0.203 sin |, It is found that (3^ = 0.434. 
Substitution into equations (22a) and (22b) then yields: 

q lfl = ijfo (0.004 - 0.0341) = -0.557 + 0.9411 

q i2 l = q 1Zo - (-0.0048 + 0.02151) = -0.0053 + 0.336i 

t 

The d amp ing factor (real part of q) in the fundamental flapping mode ' 
is here evidently negligibly affected, while the natural frequency in 
this mode is increased by about 3 percent. In the fundamental lagging 
mode here, however, the already low damping factor is dangerously 
diminished (by about-50 percent), while the already low natural frequency 
is also diminished (by about 6 percent). 

(b) Hinged blades: In this case (see Numerical example under 

BENDING VIBRATIONS IN FLAPPING and relation (17)), 

q = -0.498 + i \|(1 + 0.05 X 1.5) - (0.498) 2 = -0.498 + 0.909i 
while (see figs. 8 and 10) 

q llo = _£il^3_x_0£74 + 0> 2Y4i = -0.00623 + 0.2741 

Putting W 1 (|) = T 2 (|) = l, it is found that: ^ = O.333, 

T|. = 0.0775, anii 7l = 0.464. Moreover, assuming (|) = O.I318, 

-L L 

it- is found that = 0.0439* Substitution into equations (22a) 
and (22b) then yields: 

q lfl = ^lfo “ (°- 0016 ~ 0.01081) = -0.496 + 0.920i 


q iu - q 1Zo - (-0.00535 + 0.02541) = -0.00088 + 0.2491 


7 w 

The expressions used here for — 


approximation for 


V 
z 


(I) are based on a first 


(|) as determined by the static equation corre- 


sponding to the flapping equation (la) with the data assumed here. 
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As in example (a) for fixed "blades* the fundamental con® lex flapping 
frequency is here only slightly affected by the coupling* while in the 
fundamental lagging mode* the natural frequency and especially the 
damp ing factor are appreciably diminished. 


"VIBRATIONS OCP UNBALANCED BLADES - FLUTTER 


When the center of gravity of the cross section of a blade does 
not coincide with its aerodynamic center (Q ^ 0) then* as can be seen 
from equations (la) and (lc), the torsional vibrations will be aero- 
dynamic ally coupled with the flapping vibrations. When such coupled 
vibrations are unstable* flutter is said to occur. In the present 
section simple criterions are developed for the avoidance of the 
flutter of rotating helicopter blades. The analysis is* however* based 
on quasi-stationary flow (as in the previous sections) and is therefore* 
strictly speaking, approximately valid only for cases of low reduced 

is the magnitude of the inflow velocity. 


** w(g)e pt 

= 0(5)^* 


neglecting second-order terms there* 
it is possible to combine these equations into a single stationary 
condition. By application of the Rayleigh method* with the assumption 


frequencies ^ X ( ^ )* where "V 


5? ~ VflAR, 

Basic equations.— By putting 




0(S*t) 


* = £ 


into equations (la) and (lc), and 


*U) - «*,({) 
8(0 =»A(!) 


(23) 


for any given principal mode (characterized by subscript n), 
where I a ( i ) is the corresponding uncoupled mode of deflection in 
flapping while Y n (§) is the uncoupled mode of deflection in torsion* 
the following linear equations are obtained for the coefficients a^ 
and b n : 



30 


NACA TN 1999 




KACA TN 1999 


31 


% = MjlCI + ig 9 )f n 



H + 



^maz ~ Sste £0 r y 2 
(^max “ ^inin) 0 ° ^ 


dS - 



with 



Tlie condition for the existence of a nan-trivial solution (a^ 
and h n not "both zero) to equations (24a) and (24b) leads to the 
following quartic equation in the complex frequency q for the coupled 
mode in torsion and flapping: 

c^q 1 *- + c^q -3 + c 2 q 2 + c^ + o Q = 0 (25) 

where 

c 4 = “nVn 

°3 = “n^n + r n v n 

c 2 = J n v n + a n^n + Vn _ 

°1 = J n*n + ^ + CQ ( P n ~ °n f 
c 0 = Vn 
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For the fundamental coupled mode, the substitutions Yj_(| ) = sin 

and Xi(i) = — 44^ + 6 1 2 (fixed— ended blade) or X^(4) =• 4 (hinged 

blade) may be made for the purpose of obtaining approximate solutions 
for the c har acteristic values of q from equation (25)* 

Stability criterions .- Neglecting internal damping at first, it 
can be ft shown that the coupled torsional and flapping vibrations will be 
stable® if and only if 


Mn > 0 


(26a) 


while also 


( C 2 C 3 “ c l c k) ~ 


CqC^ > 0 


(26b) 


A simple method of taking internal damping into account after condi- 
tions (26a) and (26b) have been considered is described in appendix C. 

It may be noted that inequality (26a) corresponds to a necessary 
and sufficient condition for the stability of the uncoupled torsional 

vibrations of- unbalanced blades, that is, for the prevention of 
torsional divergence. 

It appears convenient in practice to regard inequalities (26a) 
and (26b) as conditions governing the proper- chordwise location of the 
center of gravity of a blade section with respect to the aerodynamic 
center, that is, governing the suitable values of Q or j. This is 
in accordance with the results of references 8 and 9 , where the chord- 
wise mass distribution was found to be the chief factor determining 
the possibilities of flutter of helicopter blades. 

If ,3 = 0, thatnLs, if the blades are mass-balanced, then flutter 
cannot occur, since there will be no coupling between the torsional and 
the flapping vibrations. If J < 0, then It will be found that condi- 
tions (26a) and (26b) are usually easy to satisfy, so that flutter will 
not occur. If, however, J > 0, that Is, if the center of gravity of a 
blade section is behind the aerodynamic center, then conditions (26a) 
and (26b) may be violated and flutter will occur. 

In this quasi-s tat ionary type of analysis, there are two Important 
stabilizing influences: The internal damping (proportional to gy 

and gg) in bending and in torsion and especially the aerodynamic 
damping (proportional to ^ or P) in torsion. Because of these 

8 Inequalities (26a) and (26b) are the conditions that all the real 
parts of q, as determined by equation (25), be negative, so that posi- 
tive damping occurs. 


rote 
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influences , flutter may be prevented even when the center of gravity 
is "behind, hut hot too far. behind, the aerodynamic center of a blade 
section. 

Numerical example.— Consider a fixed— ended uniform blade with the 
data in figure 2. Then, with y = £ (where € =0), 

%(!) = I 1 *" — 4^3 + 6| 2 , and Y^(£) = sin £ £ for the fundamental 

coupled mode, inequality (26a) leads to the condition J < 0.186. This 
is the condition for the prevention of torsional divergence (with 
internal damping neglected). Inequality (26b), however, leads to the 
condition 0.15 < J < 0.1 6 . Hence the critical value of j in this 
case for the prevention of flutter is between 0.15 ant. 0.l6. With 
internal damping taken into account (gg = gy = 0 . 05 ) by the method of 
appendix C, the critical value of j is found to be between 0.l6 
and 0.1T« To obtain some insight into the character of the flutter 
vibrations in an unstable case, suppose j = 0.17* Then two of the 
roots of the complete quart 1c (i.e., including the internal damping 
terms), equation (25), are found to be q = 0.123 i 3 . 971 * This pair 
of roots evidently corresponds to an unstable mode of vibration, since 
the real parts are positive, indicating "negative aerodynamic damping," 

with a logarithmic increment of - - ~ c^~^ ~ = 0.1946. The relative 

amplitudes of the flapping and torsional vibrations in this mode can 
be obtained by substituting the values found for q into equation (2ka) 
or equation (24b). It is thus found that ^ = —0.021k + 0.0076i, 

indicating that the tip amplitude of the flapping vibrations is only 

3 \j( 0.0214) 2 + (0.007 6) 2 , or 0.067 times the tip amplitude of the 
torsional vibrations. This small ratio of bending amplitudes to 
torsional amplitudes in flutter appears to be a familiar phenomenon. 

(See reference 8.) The complex value of the ratio Sq/bn illustrates, 
of course, the well-known phenomenon of phase difference between the 
coupled bending and torsional vibrations in flutter. 

Critical speed.— In addition to the design viewpoint with respect 

to chordwise mass distribution for the prevention of flutter for a 
given rotor angular speed 2, another possible viewpoint here is that 
of the critical angular speed. From this point of view, all parameters, 
including j, may be considered as given, except the rotational 
speed fl, which appears implicitly in the terms containing J n , n^, f^, 

and 7 n . Conditions (26a) and (26b) bher determine the critical range 
of 2. If j > 0, for example, then condition (26a) in actual cases 
gives an upper limit to the permissible value of fi. 
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Aerodynamic coupling between lagging and torsion . 9— It may "be 
observed that to the order of approximation used in the analysis, there 
is- no aerodynamic coupling "between the lagging and the torsional 
Titrations, since the torsion equation, equation (lc), does not contain 
any terms in the lagging displacement v, or in derivatives of v. It 
may he noted, however, that with higher-order terms such a coupling 
would exist "because of a twisting couple c^^pfirc^v, where c^ is 

the moment coefficient of a tlade section about the aerodynamic center. 
For airfoil sections for which c m ac ’f therefore, there will he 

coupling between the lagging and torsional vibrations. Since, 
however, v in equation (lc) and 9 in equation (lb) will occur as 
products with other first-order-small quantities, these terms will be 
smaller than other terms appearing in these equations, indicating 
that the coupling between lagging and torsion is in any case small. 


CONCLUSIONS 


From a theoretical analysis of the frequency and damping charac- 
teristics of the free modes of vibrations of fixed-ended and hinged 
elastic helicopter rotor blades in hovering and vertical flight, the 
following conclusions are drawn; 

1. The fundamental natural frequency in flapping of a rotating 
helicopter blade hinged at the axis of rotation is equal to the rotor 
angular speed. If the blade is fixed— ended at the axis of rotation, 
then, unless its bending stiffness is unusually high, the fundamental 
natural frequency in flapping will be only slightly higher than the 
rotor angular speed. With aerodynamic damping neglected, the rigid- 
body oscillations will be an exact mode of vibration in flapping of a 
blade hinged at the axis of rotation. 

2. The natural frequency of vibration in any principal mode in 
flapping can be expressed approximately by means of a simple equation. 
In this equation the centrifugal contribution to the square of the 
natural frequency of a rotating blade in any mode is virtually the 
same whether the blade is hinged or fixed at the root. However, the 
elastic contribution depends on the boundary conditions at the root. 

A simple method of calculating the natural frequencies of a rotating 
beam in bending, based on these considerations, was developed. 

3. A sinple relation exists between the lagging and the flapping 
natural frequencies of any undamped uncoupled principal mode of 
vibration. 

%he discussion in this paragraph is valid whether the blades are 
mass— balanced or not. 
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4. If the "blades are attached at a ‘small distance (eccentricity) 
away from the rotor axis of rotation, then e.n the natural flapping and 
lagging frequencies are increased. This effect is especi ally important 
in the lagging motion of hinged "blades, where there would otherwise "be 
no restoring forces in the fundamental mode. 

5 . The flapping vibrations are heavily damped aer odynami cally in 
the lower modes, especially the fundamental. The internal damping in 
these modes is, in fact, negligible in comparison with the aerodynamic 
damping. 

r 

6 . The aerodynamic damping in lagging is proportional to the 
profile-drag coefficient of a blade section and is much smaller than 
that in flapping. Consequently, the internal damping in the lower 
modes of vibration in lagging is as important as the aerodynamic 
damping. 


7. The aerodynamic logarithmic decrements in flapping, lagging, 
and torsion diminish with the principal modes, in approximately inverse 
proportion to the vibration frequency. 

8 . In flapping and torsion, the importance of aerodynamic damping 
relative to that of internal damping diminishes with the principal 
mode of vibration. In lagging, however, this is not quite true in the 
lower modes, since for fixed-ended blades of relatively low bending 
stiffness in lagging the internal logarithmic decrement decreases 
sharply from the fundamental to the second mode. 

9 . The centrifugal torque In a rotating helicopter blade usually 
exerts a negligible effect on the natural torsional frequencies. There 
is considerable aerodynamic damping in torsion, due to the rotation of 
a blade. 

10. The effect of the free torsional vibrations on the flapping and 
lagging vibrations of mass-balanced blades is in practice negligible, 
because of the relatively high natural frequencies in torsion. 

11. In general, Coriolis, aerodynamic, and centrifugal coupling 
forces exist between the flapping and lagging motions of a rotating 
helicopter blade. These forces have a second— order effect on the 
flapping and lagging natural frequencies and damping factors, which is 
nevertheless appreciable in the fundamental lagging mode. 
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12. If the blades ore unbalanced, then the resulting coupled 
flapping and torsional vibrations may be unstable if the cross-sectional 
center of gravity is too far behind the aerodynamic center. Simple 
criterions for prevention of such a type of flutter, based on quasi— 
stationary flow, were derived. 


Polytechnic Institute of Brooklyn 
Brooklyn, N. Y., May 2, I 9 I 4.9 
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APPENDIX A 

TORSIONAL VIBRATIONS OF NCJTUNI3PORM BLADES 


Equation ( 5 ) can be mathematically transformed into the following 
stationary condition: 



m 4 8 ' 


+ Mg 


^max 


- I. 


■min 


^rrdn) { 


J o 


d| - 0 


(Al) 


where 5 is here an operator denoting variation. By assuming any mode 
of deflection to he expressible in the form of the series 


8(1)- H !>* sin k § I • (82) 

k— 1 j3j5j * • • 


the following set of linear homogeneous equations* written in matrix 
form, is obtained in the coefficients b^-: 



with m = 1 , n = 1 , and 
constants defined thus; 


“W = ^max* 801 where ^ and v mn are 
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T im = M 1 I ^ ^ — cos H 4 cos ^ 4 d| + 


Pi 


„ I W ^min mit „ nir „ 

Mg | — r — sin — 4 sin — 4 d| 


(l-max ^min) c 


v 


mn 



. nut „ . 

sin — | sin 
2 


nit 

2 


t d| 


It will "be found in actual cases that the 1*^— term in the expression 
for is negligible in comparison with the Mj— term. 

Equations (A 3 ), written out explicitly, look thus: 

^l^ll + i 2v n) + ^(^l + ^ V 3l) + * * * + + &Bl) + • • • = 0 

+ ^ 13 ) + * 3 (*33 + ^ V 33) + • • • + + ^jqp) + • • . = 0 

^rn + ^In) + + 1 2 v 3 n) + * * * + *n(V + ^nn) + • • • - 0 

The condition for the existence of a non— trivial solution to 

equations (A3) is that the determinant of the coefficients of- these 
equations vanish. This leads to the following determinant equation: 


-^nm + 1 v mn 


= 0 


(M) 


Equation (A4-) determines the natural frequencies of the various 
principal modes. Although equation (At) is, for exact calculations, of 
infinite degree with an infinite number of terms, it will he found in 
actual calculations that this equation is rapidly convergent; that is, 
that it suffices to use only a few terms in series (A2) to obtain to a 
sufficient- approximation the values of the natural frequencies of the 
lower modes. This convergent property is one of the chief advantages of- 
the Bayleigh-Ritz method- outlined here. 
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APPENDIX B 


MATHEMATICAL DETAILS FOR BENDING- VIBRATIONS IN FLAPPING- 

(a) Proof that relation (11) would "be exact if and only if the 
centrifugal loads -would hare no effect on the modes of deflection of a 
rotating "beam: 

Suppose , first, that relation (11) is an exact general relation; 
that is, suppose that the relation 


? 9 ? 

In = len + len 


(Bl) 


is in general exactly valid. By definition of and q Qn it 

follows fro.'i equation (10) that 




Aq cn 








(B2) 


where w cn (|) is a principal mode (characterized hy subscript n) of 
deflection of a rotating cable without bending stiffness, while w en ( £) 

is the corresponding mode of deflection of a nonrotating beam. By 
putting equations (Bl) and (B2) into equation (10), the following 
equation is obtained: 



v n Al 
v en \Iio 



^( TW cn , )‘ ( Tw n* )* 


= 0 


(B3) 


where w n ( | ) is the mode of deflection of a rotating beam with bending 
stiffness. 
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In order that relation (Bl) he exact regardless of the value 
of Ki it is necessary that equation (B3) he satisfied for any value 
of and since w en (|) and w cn (|) are independent of Kq, it is 

therefore necessary thatrthe following two equations he satisfied: 



- M - 0 (B5) 

The unique solution of equation (B4) for w^i- ) satisfying all 
of the boundary conditions is 

w n^ )=w en^) ( b6 ) 

Similarly the solution of equation (B5) for v n (| ) Is 

' w n^ =w cn^ ^ B7) 

Equation (B 6) proves that- relation (Bl) could he exact only if the modes 

of deflection of a rotating "beam were exactly the Bame as those of the 
same "beam when it is stationary. 

To prove the converse, suppose relation (b 6) were In general 
exactly valid. Then ‘substitution of relation (b 6), together with equa- 
tion (B2), into equation (10) would lead to the following relation: 

^ (<ln 2 " W^en^) " ‘W^cnU)] + [(™cn*)' " ( TW en*)*] = 0 ( b8 ) 

Observing that w en (t ) and- w cn (| ) are functions not containing q Qn 2 
i_t follows from equation (B8) that w cn (^ ) = w en (| ), while 

also q n 2 = + len 2 ' Tilus relation (b 6) would Imply relation (11). 
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(b) Proof that relation (11) cannot he an exact general statement, 
vh.eth.er the hlad.es are hinged or fixed— ended: 

It has been shown that relation (11) implies relation (B7). It 
will he shown, however, that relation (B7) cannot he an exact general 
statement because w^d) cannot in all cases satisfy the boundary 
conditions v cn " = = 0 required at the free tip of a beam with 

bending stiffness. 

By definition, v on (5) satisfies the equation (cf. equation (10)): 
-( TW cn')‘ + a -cn 2w cn = 0 ( B 9) 


Differentiating equation (B9) three times with respect to |, noting 
that t* = — | (for simplicity, it is assumed here that ^ , and 


assuming that = w cn m =0 at 5=1, the following relation is 

obtained, where vq, and derivatives thereof denote the values 

of v cn (5), of (|), and of their derivatives at 5=1: 


■Ao 




w^’ + 



(BIO) 


Since w cn (5) = w en (5), it follows from equation (B2) that 




(Bll) 


Moreover, equation (B9) implies 

v l* = “^cn^l 


(B12) 


By putting equations (Bll) and (B12) into equation (BIO), the following 
relation, required to be valid for any type of blade in every principal 
mode, is obtained: 
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*G8,V ■ 3 (^ 2 * (B13) 

It is evident that relation (B13) cannot "be valid in general 10 since, 
for example, in the case ( —) ^ 0 equation (B13) leads to an 

Wi 

absurdity, the left side being negative while the right side is positive 
or zero. 

(c) Proof that there are no exact values of q^ 2 for a fixed— 
ended blade: 

If there were exact values of q^ 2 then these would satisfy 
equation (B9) , with the appropriate functions v cn (|). However, for 
fixed— ended blades, v cn (0) = w CI1 , (0) = 0. Hence from equation (B9) 
it follows that w cn (£ ) would have to satisfy also the condition 
tion w cn "(0) = 0. Moreover, by successive differentiation of equa- 
tion (B9), it is seen that all of the derivatives of w cn (| ) would 
have to vanish at the root of the blade. Hence w cn (| ) = 0. This 
indicates that there cannot be any exact values of q^ 2 with corre- 
sponding functions w cn {£ ) satisfying equation (B9) for a fixed-ended 
blade . 

(d) Numerical check on accuracy of approximation of relation (11): 

Equation (10) can be mathematically transformed into the following 
stationary condition: 

5 J | K i ^ v ” 2 + | TV ’ 2 + 1 ^ i 2 * 2 d i = 0 cb ^) 

The Eayleigh-Pitz method can be applied to condition (B14) by setting 

w= 28^(6) (B15) 

where 2^(1) are given functions satisfying the boundary conditions. 

10 There may be special cases in which equation (BI3) is valid. 

For example equation (B13) may, for a given blade, be valid in one 
particular mode, as in the fundamental mode of a hinged blade, 
where q e q 2 = 0 while - — 1* 
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Tiie most general sets of polynominals independently satisfying tits 
"boundary conditions for fixed— ended and for hinged "blades are: 


fixed-ended: 


| m - i m(m - l)(m - 2)£ 3 + i m(m - l)(m - 3)£ £ 
m ^ k 


hinged: 


Z^Ci) = i m + | m(m - l)(m - 4)| 3 - jg m(m - l)(m - 3)1* 


y (bi6) 


11 “ 5j Tj • • • 


The Rayleigh-Ritz method (see* for example , reference 10) was 
applied here for fixed-ended blades with the first three terms of equa- 
tions (Bl6) in series (B15). By substituting the series for w into 
the integrand of condition (Bl4) and differentiating this integrand 
with respect to each of the three coefficients a^, a set of three linear 
homogeneous equations in a^ was obtained, the vanishing of whose 
determinant led to a cubic equation in the negative squared natural 
frequencies q 2 . By investigating the solutions of this cubic, it was 
found that all the three roots could, to a high degree of approximation, 
be expressed in the form q^ 2 = -fgnEq - f cn where f^ and f cn were 
positive numbers, independent of Kq. This verified relation (11), 
since the quantity -fgn^l could then be interpreted as q^ 2 , while 

the quantity - f cn could be interpreted as q^ 2 . 


(e) Effect of root conditions on centrifugal contributions to 
natural frequencies : 

As can be seen from reference 5 (appendix), the values of q cn 2 
in relation (11) are independent of whether a blade is hinged or fixed 
at the root, since the only boundary condition which can in an exact 
solution, by Bessel series, be satisfied by equation (B9) for w^d) 
is w cn (0) = 0. It is significant to check whether the approximate 
methoa outlined in item (d) leads to the same results in actual calcu- 
lations. For this purpose the stationary condition, equation (Bl4), 
with Kq = 0 was applied by using four terms in series (B15) with X^d ) 
for fixed-ended blades and then with Z^l) for hinged blades (see 
equations (Bl6)). For fixed-ended blades the following values of q^ 2 
were obtained for the first two principal modes; q<,q 2 = —1.035 

o 

and 


<*c2 


2 _ _i 

. IlOU. x WJ. i J.J. w w WH V t '3-CX 

2 = —6.21. For hinged blades, the corresponding values obtained 
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were q^ 2 = — i 1c2 2 = -6.05- Thus the values of q cn (n = 1,2) 

obtained by using series (Bl6) for fixed-ended blades were within 
2 percent of those obtained by using the series for hinged blades. 

As a further numerical check the Rayleighr-Ritz method was applied 
to condition (Bl4) with Kp = 0 with the first four terms in 
series (B15) with 


X m U) 




(BIT) 


This series satisfies only the single boundary condition of zero 
deflection at the root. The roots obtained for q^ in the first 

two modes were: q Q ^ 2 = —1 and q^ 2 = -6. Thus the simple series, 

equation (BIT ) , led to nearly the same results as either of series (B 15 ). 
It can, in fact, be shown that for a constant or linearly varying cross- 
sectional area the values of q^^ thus obtained in any mode will be 
the same as the exact values obtained by Bessel series, provided a 
sufficient number of terms is used in series (BIT) • 

(f ) Method of calculating natural frequencies of any mode for 
rotating blade of variable cross section: 

By making use of relation (11) and of the fact that the centrifugal 
contributions q cn 2 to the negative squared natural frequencies are 
negligibly affected by the root conditions of- a blade, the Raylelgh-Eltz 
method applied to condition (Bl4) leads to the following approximate 
method of calculating the elastic, negative squared frequencies q^ 2 . 

The values of q en 2 are the roots for q 2 of the determinantal 

equation: 11 


11 The determinant in equation (Bl8) is obviously diagonally 
symmetric. 



NACA TE 1 999 


^5 


JllMllI 2 

J l5' Kt l5 ( l 2 

J n 6^6^ 

2 

* ■ * J lk +a lk <1 

Vta 5 - 2 


j 56"56^ 

* • * J 5k +a 5k < ^ 2 

J eL«a.^ 


J 66 +a 66 q2 

* * * J 6k w 6k!^ 

' J kl +< \l (l2 



Jkk 4 ^ 9 - 2 


where 




(|)Z ll ,, (|)X i "(|) dl 


(Bi8) 


“hi- 



A (Dx^OXiCi) as 


and or X^(§) is given By series (Bl6), which satisfies all of 

the Boundary conditions in Bending. -^- 2 

The values of q^ 2 are the roots for q 2 of the determinantal 
equation: 


T ii K1 ii < i 2 

T 12 +a l2 q2 

Il 3 «l 3 l 2 

. . . T iz +a n < l 2 

^l^l*! 2 

1'22 +C£ 22^ 2 

T 23 4<i 23 q 2 

• * ’ T 2l +a 2z q2 

T 3i*°3i^ 

i ^32"*^ :£ 32 1 ^ 2 

T 33 m 33 !l2 

. . . ^3 l +a 3l^ 2 

T 21 +a 21« 2 

T l2 +a Z2 <l2 

T J3 +a ^3 9 - 2 

. . . T jz + “n q2 


12 If the Blades are fixed-ended instead of hinged at their roots, 
then, according to series (Bl6), the suBscript 1 in the determinant in 
equation (Bl8) should be replaced By the suBscript h. 
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where 


with, 


T f g $ 


T(6)x f '(l)x, 


g 


(«) <14 


t(i) S 



<e> T « 


and X f (|) or X g (|) is now given "by the simple series, equation (B17) 
which satisfies only the condition of zero deflection at the root. In 
equation (Bl8), (k — 3) is the number of terms used in series (Bl6), 
while in equation (B19), l is the number of terms used in series (B17) 


The chief advantage of the method, outlined is that the convergence 
is rapid, even for modes above the fundamental. Consequently, only 
relatively few terms need be used in actual cases in series (B15) to 
obtain the natural frequencies to a sufficient approximation. Suppose, 
for example, it is desired to obtain the natural frequencies for the 
three lowest modes. Then one first chooses k = 3 and l = 3. The 
three values of q 2 obtained from equations (Bl8) and (B19) represent 
the values corresponding to the three lowest modes. The two lowest 
values of q 2 will usually be more accurate than the third. To improve 
the accuracy one then chooses k = 4 and 1=4 and obtains again the 
roots for q 2 corresponding to the first three modes, as well as an 
additional root corresponding to the fourth mode. It will usually be 
found that at least the two lowest roots for q 2 thus obtained are 
negligibly different from the two lowest roots obtained with k = 3 
and l = 3. This means that the natural frequencies of at leash the 
two lower modes have been determined with sufficient accuracy. The 
natural frequency of the third mode could be obtained to a still further 
approximation by taking k = 5 aad l =_5 • In general it can be stated 
that, if for k = k (say) and for k = k + 1 the s lowest roots 
for q 2 of equation (Bl8) are practically the same, then these are the 

values of "p — 

**1 

valid for q^ 2 


for the s lowest modes. The method is similarly 


If the order of the determinant in either equation (Bl8) or (B19) 
is higher than three, then experience has indicated that a convenient 
method of solving the determinantal equation for any root is by 
systematic trial and error. One chooses two values of q 2 which are 
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telieved to "be fairly close to an actual root and evaluates tie deter- 
minant for each of these values. The value of q 2 then obtained by 
linear interpolation (preferably) or extrapolation, according to which 
the determinant would be zero, will then be a closer approximation to 
the actual root. 

A convenient method of evaluating a deter minant of any order is 
to transform it into a triangular form. This method will be illustrated 
here for a third— order determinant. Consider the determinant 


a ll 

*21 

a 31 

a „ 

a 

a 

12 

22 

32 

a 13 

"23 

a 33 


It is desired to obtain a determinant equal to D where the elements 
in the present position of a^, a^ 1 , and a^ 2 will be zero. 

Multiply the first column by a 2l /a 11 and then subtract this column 


from the second column. Also, multiply the first column by a^ /a, . 

and subtract it from the third. Then a determinant of the following 
form is obtained: 


a ll 

0 

0 

a i2 

CVJ 

CM 

& 

*32 

a 13 

*23 

*33 


In D* now multiply the second column by bg 2 /t>22 81111 subtract It from 

the third column. Then the following result is obtained: 
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Th.es 0 calculations should he performed with a sufficient number of 
significant figures to avoid any large errors due to relatively small 
differences . 
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AEPSNDH C 


INTERNAL DAMPING IN UNBALANCED BLADES 


The characteristic equation* equation (25)* including internal 
damping can he "written in the form: 

F 1 (q) + iE 2 (q) = 0 (Cl) 

where Fq and F2 are polynominals in q with real coefficients. 

The effect of internal damping is represented by F2. 

Since |F 2 (q)| « (N’ 1 (q.) J ^ equation (Cl) can he solved for q to a 
good approximation by applying Newton’s method as follows. The equation 

*q(<l) = 0 (C2) 

can he solved first. Suppose one of the roots thus obtained is qQ. 

Then a corrected value qq of q 0 due to the internal damping will he 


iF 2 (q Q ) 

q x = q Q (C 3 ) 

Fi’(qo) + O’ 2 , ( ( lo) 

Equation (C3) can he used to determine the critical value of j 
thus. Suppose that from condition (26h) it is found that to prevent 
flutter it is necessary that j < (say). Then J Q is the critical 
value of j with internal damping neglected. To take internal damping 
into account, slightly higher values of J may he chosen, and the 
corresponding roots of the quartic, equation (25), including internal 
damping can then he obtained hy means of equations (C2) and (C3). The 
minimum value of J for which at least one root of the complete quartic, 
equation (25), will have a positive real part will he the critical 
value of j with internal damping. 

Near the critical value of 3, the imaginary part of q Q will 
generally he much greater than the real part, so that the actual " 
computations can he simplified hy putting q Q = icu (ca real) in 
equation (C3). 
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Steel spar 



'-Wooden rib web (a m A = - 4 - ,) 


n = 25 radians/sec - 239 rpm 
R = 20 ft 
W = 3000 lb 

mi- i. t _ j . _ 

jLnree uiaues 

d/c = 1/8 


Figure 2.- Blade section used in numerical examples 
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Figure 5.- Aerodynamic and internal logarithmic decrements of hinged blade 
in flapping. Constant cross section; K-j_ = 0.00400; e = 0.05. 
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Figure 6.- Aerodynamic and internal logarithmic-decrements of fixed-ended 
blade in flapping. Constant cross section; = 0.00400; e = 0.05. 
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and Kg 

Figure 7.- Variation of co en 2 with K^_ (flapping) or K 2 (lagging) for any 

mode (n). 
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Figure 9.- Lagging frequencies of fixed-ended blade. Constant cross section; 

Kg - 0.00400; e = 0.05. 
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Principal mode. of vibration 


Figure 10.- Aerodynamic and internal logarithmic decrements of hinged blade 
in lagging. Constant cross section; Kg = 0.00400; e = 0.05. 








